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^  Dispersion4"  effects  are  considered  in  addition  to  "Location44 
effects  of  factors  in  the  inferential  procedure  of  sequential 
factor  screening  experiments  with  m  factors  each  at  two  levels 
under  search  linear  models.  Search  designs  in  measuring 
,  '^Dispersion""  and  "Location*  effects  of  factors  are  presented  for 
both  stage  one  and  stage  two  of  factor  screening  experiments  with 


4  <  m  <  10. 


1.  INTRODUCTION 

In  a  factor  screening  experiment  the  problem  is  to  find  a  few 
effective  factors  from  a  list  of  a  large  number  of  possible 
factors  influencing  the  response  and  screen  out  all  non-effective 
factors.  Watson  (1961),  Patel  (1962)  and  others  considered  only 
"location44  effects  in  their  study  of  influence  of  factors  on 
response.  Taguchi  and  Wu  (1985),  Box  and  Meyer  (1986)  and  others 


■r-  J  ^  -  ■  -**-  s  ■  -*■’*.*■  ^  V  *■*  * ■r-  *-" -  ■.*-*>  V- v—  -r .  J t.  -*  A “>  J  A  ^  ..  J 


86  12  1 1 


OG 


considered  both  "location"  and  "dispersion"  effects  in  the  study 
of  influence  of  factors  on  response.  The  concept  of  "dispersion" 
effects  in  factor  screening  experiments  was  introduced  in  Taguchi 
and  Wu  (1985)  for  replicated  fractional  factorial  designs  and  in 
Box  and  Meyers  (1986)  for  unreplicated  fractional  factorial 
designs.  Srivastava  (1975,  1976)  introduced  the  search  linear 
model  and  showed  that  factor  screening  problems  can  be  solved 
using  designs  under  the  search  linear  model  (called  search 
designs)  which  have  smaller  numbers  of  treatments  than  that  in 
comparable  designs  in  the  literature.  Ghosh  (1979),  Ghosh  and 
Avila  (1985)  presented  many  such  search  designs.  However,  those 
search  designs  were  constructed  to  measure  the  "location”  effects 
only  and  most  of  them  are  unable  to  measure  the  "dispersion" 
effects  of  factors.  In  this  paper  we  define  the  "dispersion" 
effects  of  factors  under  the  search  linear  model  and  present 
search  designs  to  measure  both  the  "location”  and  "dispersion" 
effects  for  factor  screening  experiments. 

Throughout  this  paper,  we  consider  the  sequential  factor 
screening  experiments  with  n  factors  each  at  two  levels  under 
search  linear  models.  At  the  first  stage  of  the  experimentation 
the  problem  is  to  estimate  the  dispersion  effects  at  two  levels  of 
factors  under  search  linear  models,  to  use  the  dispersion  effects 
in  finding  the  most  effective  factor  out  of  m  factors  and  finally 
to  determine  the  optimum  level  combination  of  the  most  effective 
factor  using  the  "signal  to  noise"  ratio.  At  the  second  stage, 
the  problem  is  to  estimate  again  the  dispersion  effects  at  four 
level  combinations  of  everv  two  factors  under  search  linear 
models,  to  find  two  most  effective  factors  out  of  m  factors  in 
presence  of  interactions  anc  finally  to  determine  the  optimum 
level  combinations  of  two  most  effective  factors.  The  process 
continues  until  we  find  all  effective  factors.  In  section  2  of 
this  paper  we  discuss  the  models  and  inferential  procedures.  In 
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section  3  we  present  search  designs  (A  £  m  ^  10)  for  both  stage 
one  and  stage  two  of  factor  screening  experimentation. 


2.  INFERENCE 

In  a  2ni  factorial  experiment,  the  treatments  are  denoted  by 

(a  ,,..., a  ),  a. =0,1;  the  general  mean,  main  effects  and  two  factor 

interactions  are  denoted  by  p,  and  F^F^,  respectively;  the 

observation  corresponding  to  the  treatment  (a  , . . .  ,a  )  is  denoted 

1  m 

by  y(a  ,...,a  )•  The  expectation  form  of  the  model  is 
I  m 


E(y(a  . a  ))  =  p  +  £  b  F  +  Z  I  b.  b.F,  F.+. . . ,  (1) 

1  m  1  i=l  i  i  i=l  j=l  1  j  1  j 

i<j 


where  =1  if  a^=l  and  b^=  -1  if  a^=0. 
Stage  1 


At  the  stage  1  of  the  experimentation  we  consider  m  different 
models  and  the  ith  model  is 

y (a  ^  .  ,a  ^ .  ,a^) )  =  y  +  b_^F^,  i=l,...,m,  (2) 

V(y(a1,...,ai,...,am))  =  (3) 

and  the  observations  are  uncorrelated,  where  the  variance  o^(0) 

2 

and  o  ^1),  i=l,...,m,  are  unknown  constants,  called  the 

"dispersion"  effects  for  the  factor  i.  We  now  consider  a 

fractional  factorial  design  with  Nj  treatments.  (Treatments  may 

or  may  not  be  all  distinct.)  We  denote  for  u  =  0,  1, 
u  u 

y^  =  the  vector  of  observations  corresponding  to  the  treatments 
with  the  ith  factor  at  the  level  u, 

y.  =  the  simple  arithmetic  mean  of  FU  observations  in  y?, 

1  i.  — 1 


s|(u) 


K  - 1) 
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Note  that  N. 


2  (N°  -  l)  sj(0)  +  (N*  -  l)  Sj(l) 

‘‘  (Nj  -  n)  +  (N‘  -  1) 

N?  +  N* .  We  now  write  (2)  and  (3)  as 
l  i 


1  N?  '  JN? 


1  ,  • 


b  -  F 


W  +  F 


oj(0)l 


W  V  *V 

where  j  is  a  vector  with  all  elements  unity.  It  can  be  checked 
that  E(s^(u)J  =  o*(u)  and  the  generalized  least  squares  (and  also 
the  ordinary  least  squares)  estimators  of  (b-F^)  and  (b+F^)  are 

y°  and  y*,  respectively.  We  select  the  ith  factor  as  the  most 
effective  if  s2  is  a  minimum  for  i  e{l,...m}.  We  select  the  uth 
level  of  the  factor  i  as  the  optimum  level  if  the  (square  of) 
"signal  to  noise"  ratio  [y?/si(u)J2  is  a  maximum. 

Stage  2 

At  the  stage  2  of  the  experimentation  we  consider 
different  models  and  the  (i,j)th  model  is 

E(y(a1,...,ai,...,aj,...,am))  =  b+b^F^+b  ^F  j+bj,b^.F^F  ^ , 

i,j=l,...,m,i^j,  (6) 

v( y( (a j » • • • > • • • >a j > • • • »an ) )  =  (?) 

and  the  observations  are  uncorrelated,  where  the  variances 
o^j(OO),  o2j(01),  o2j(10)  and  a2^(ll)  are  unknown  constants, 
called  the  "dispersion"  effects  for  the  factors  i  and  j.  We  now 
consider  a  fractional  factorial  design  with  N2  treatments.  Note 
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again  that  a  particular  treatment  may  or  may  not  be  replicated  in 
the  design.  We  denote  for  u,  v=0,l, 

UV  yy 

— ij  *  the  vector  observations  corresponding  to  the  treat¬ 

ments  with  the  levels  of  factors  i  and  j  are  u  and  v, 
respectively , 

simple  arithmetic  mean  of  N. observations  in  y?V» 


s1j(uv)  = 


UV  UV  .  \  '  /  U V  UV  . 

13  slL 

KI  -  0 


1  1 


Z  Z  (N^  -  1)  b2  (uv) 

2  u=0  v=0  J 

s±.  -  -  _ 

Z  Z  (N^  -  1) 
u=0  v=0  J 

Observe  that  N2  =  N°°  +  +  n“.  We  write  (6)  and  (7)  as 


ij 
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The  orthogonal  fractional  factorial  designs  are  indeed 
efficient  but  have  restriction  on  the  number  of  treatments  in  the 
design.  The  designs  used  in  Taguchi  and  Wu  (1985),  Box  and  Meyer 
(1986)  are  all  classical  orthogonal  designs.  In  this  section  we 
construct  designs  to  measure  the  "dispersion"  effects  defined  in 
section  2  for  the  stage  1  and  stage  2  of  the  experiment,  relaxing 
the  restriction  on  the  number  of  treatments.  We  want  to  have 
enough  observations  (at  least  two!)  in  measuring  the  "dispersion" 
effects,  the  number  of  treatments  to  be  small  and  furthermore  the 
designs  to  be  near  orthogonal. 
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3.1  DESIGNS  FOR  THE  STAGE  I  EXPERIMENT 

The  design  condition  obtained  in  Srivastava  (1975)  requires 
that  for  all  i  and  j,  i*j,  i,j  in  {l,...,m},  p,  and  F  are 
unbiasedly  estimable  (u.e.)  under  the  model  (1)  assuming  all  other 
parameters  to  be  zero.  We  denote  a  design  by  a  treatment  matrix 
T(N|Xm)  with  rows  as  treatments  and  columns  as  factors. 

Theorem  1.  A  necessary  and  sufficient  condition  that  a  treatment 
matrix  T(Nixm)  satisfies  the  design  condition  is  that  for  every 
submatrix  T1(N1x2)  of  T,  at  least  three  pairs  out  of  the  four 
pairs  (00),  (01),  (10)  and  (11)  appear  as  rows  in  TA. 

Proof.  To  estimate  the  three  parameters  p,  F^  and  F^,  we  need 
three  independent  equations  in  parameters  from  (1)  under  the 
assumption  that  all  other  parameters  are  exactly  equal  to  zero. 

For  the  submatrix  Tx(Nxx2)  of  T  corresponding  to  the  factors  i  and 
j,  any  three  rows  out  of  four  possible  distinct  rows  (00),  (01), 
(10)  and  (11)  will  give  three  independent  equations  in  p,  Fj  and 
F  under  (1).  This  completes  the  proof. 

The  characterization  in  Theorem  1  is  so  simple  that  the  checking 
can  even  be  done  by  eye  inspection.  A  result  similar  to  Theorem  3 
is  also  available  in  Srivastava  (1975). 

We  now  present  some  designs  for  4  m  _<  10  satisfying  the 
condition  in  Theorem  1.  In  the  Table  I,  for  brevity,  we  indicate 
a  treatment  (i.e.,  a  row)  in  T  by  the  positions  where  the  level  1 
is  occuring.  The  treatment  (0,...,0)  will  be  denoted  by  0.  To 
illustrate  this  for  m=4,  T( 5x4)  matrix  is  given  below. 

1111 
0 
1 
0 


T  = 


n 

1 

0 

1 


0 


0 

0 

1 

1 


0 

0 

1 

0 


We  denote  this  T  in  our  representation  as  (1234,  0,  12,  34,  13). 
Notice  that  in  every  column  of  T,  u(=0,l)  is  appearing  either 
twice  or  thrice. 


TABLE  I. 

DESIGNS  (4  <  m  <  10)  FOR  THE  STAGE  1  EXPERIMENT 


The  design  condition  obtained  In  Srivastava  (1975)  requires 

for  all  i,  j  (i*j),  u,  v  (u*v),  (i,j)  *  (u,v),  i,  j,  u,  v  in 

fl,...,m},  u,  F,,  F.,  F.F.,  F  ,  F  and  F  F  are  u.e.  under  the 
1  i  y  i  j’  u’  v  u  V 

model  (1)  assuming  all  other  parameters  to  be  zero. 

Theorem  2.  Suppose  a  design  ^^xm)  is  such  that  p,  ,  F^ ,  Fy 

are  u.e.  for  every  distinct  i,  j  and  u  in  {l,..,m}  under  (1) 

assuming  all  other  parameters  except  F.F.  and  F.F  to  be  zero. 

i  3  i  u 

Then  F.F  and  F.F  are  also  u.e. 
l  j  i  u 

Proof.  It  can  be  seen  using  Theorem  1  that  F.F  +F , F  and 

i  j  l  u 

-F  F  +F  F  are  u.e.  This  completes  the  proof, 
i  j  i  u 

Theorem  3.  Suppose  a  design  T  (N2Xm)  is  such  that  p,  F^,  F^,  F^ 
and  F^  are  u.e.  for  every  distinct  i,  j,  u  and  v  in  {l,...,m} 


under  (1)  assuming  all  other  parameters  except  F  F  and  F  F  to  be 

i  j  u  v 

zero.  Then  F  F  and  F  F  are  u.e.  if  and  only  if  for  Ti(N2x2) 
i  j  u  v 

corresponding  to  factors  (i,j)  and  T2  (N2X2)  corresponding  to 
factors  (u,v)  there  are  two  distinct  rows  in  T  with  the  number  of 
l's  in  and  T2  as  follows. 


Row 

The  number  c 

if  l's  in 

T 1 

t2 

1  a 

even 

even 

or 

b 

odd 

odd 

2  a 

odd 

even 

or 

b 

even 

odd 

Proof.  It  can  be  seen  that  from  the  rows  1  and  2  that  F  F.+F  F 

i  J  u  v 

and  -F.  F.+F  F  are  u.e.  This  completes  the  proof, 
i  j  u  v 

3.2.1  CONSTRUCTION  OF  DESIGNS  m=4 

The  treatment  matrix  consists  of  9  treatments  as  rows;  the 

first  eight  treatments  are  solutions  of  xi+X2+X3+xi4=0  over  the 

finite  field  GF(2)  and  the  last  treatment  is  (1000).  The 

treatment  matrix  is  represented  as  (1234,0,12,34,13,24,14,23,1). 

The  first  eight  treatments  form  an  orthogonal  resolution  IV 

plan.  It  follows  from  Theorem  2  that  p,  F,  ,  F.,  F.  F.,  F  ,  F  F 

i  j  ij  u  iu 

are  u.e.  for  any  three  distinct  i,j  and  u  in  {l,...,m}.  It  can  be 

checked  that  p,  F_,  F.,  F  ,  F  ,  F.F.+F  F  are  u.e.  for  any  four 
i  y  u’  V*  i  j  u  V 

distinct  integers  i,j,  u  and  v.  It  can  be  seen  that  for  all  rows 

in  the  treatment  matrix  the  numbers  of  l's  in  submatrices 

corresponding  to  (i,j)  and  (u,v)  are  either  (even,  even)  or  (odd, 

odd)  and  therefore  the  treatment  matrix  does  not  satisfy  the 

conditions  in  Theorem  3.  When  the  treatment  (1000x  is  added  to 

the  treatment  matrix,  it  follows  from  Theorem  3  that  p,  F^ ,  F  , 

F  F.,  F  ,  F  and  F  F  are  u.e. 
i  j  u  v  u  v 


The  treatment  matrix  consists  of  11  treatments  as  all  (1x5) 
vectors  with  two  elements  are  unity  and  the  other  elements  are 
zero  and  furthermore  the  treatment  with  the  levels  of  all  factors 
are  unity. 

This  is  a  balanced  array  of  full  strength.  There  is  a 

complete  2  factorial  experiment  w.r.t.  any  three  factors  and 

4 

there  is  a  resolution  V  plan  for  a  2  factorial  experiment  w.r.t. 
any  four  factors.  Thus  the  design  satisfies  the  design  condition. 
m=6 

Design  1:  The  treatment  matrix  consists  of  16  treatments  as 
follows : 

The  14  treatments  are  solutions  of  xi+X2+X3+X4=0  and 
X3+X4+X5+xg=0  over  the  finite  field  GF(2)  excluding  the 
treatments  (000000)  and  (111111);  and  the  other  two  treatments  are 
(10  00  00)  and  (00  00  01). 

3 

It  can  be  seen  that  all  2  distinct  treatments  are  present  in 

rows  of  T  w.r.t.  any  three  factors  F.,  F.  and  F  .  It  thus  follows 

i  J  u 

p,  F.,  F.,  F.F.,  F  and  F.F  are  u.e.  for  every  three  distinct 
’  i  j  i  j  u  l  u 

factors  i,  j  and  u.  For  four  factors  of  the  type  (F| ,F2  ,F3  ,F4 ) , 
(F| ,F2 ,Fg ,Fg )  and  (F3 ,F4 ,F5 ,Fg) ,  we  get  basically  the  design  for 

Lf 

m=4  with  replications.  We  get  a  resolution  V  plan  for  a  2 
factorial  experiment  w.r.t.  any  four  factors  other  than 

(Fi ,F2 >f3 »f4 ).  (Fi >f2»f5»f6)  and  (F3 »F4 >F5 >F6 ) •  Therefore, 

p,  F  ,  F  ,  F  F  ,  F  ,  F  and  F  F  are  u.e.  for  every  four  distinct 
ljijuv  uv 

factors  i,  j,  u  and  v. 

Design  2 :  The  treatment  matrix  consists  of  15  treatments  as 
follows . 

The  15  treatments  are  solutions  of  x i+X2+x3+x4=0  and 
x3+x4+x5+x6=0  over  the  finite  field  GF(2)  excluding  the  treatments 
(01  01  01),  (10  10  10)  and  (00  00  00);  and  the  other  two 
treatments  are  (10  00  00)  and  (00  00  01). 


The  argument  is  similar  to  that  in  Design  1 . 

m=7 

The  treatment  matrix  consists  of  17  treatments  as  follows. 

The  14  treatments  are  solutions  of  x A+X2+X3+X4=0, 

X3+xlt+X5-tocg=0  and  X|+X3+X5+X7=0  over  the  finite  field  CF(2) 
excluding  the  treatments  (1111111)  and  (0000000);  and  the  other  3 
treatments  are  (10  00  00  0),  (00  00  01  0),  (00  00  00  1).  [One  may 
keep  (1111111)  in  the  design  although  the  design  is  all  right 
without  it . J 
m=8 

The  treatment  matrix  consists  of  18  treatments  as  follows. 

The  14  treatments  are  solutions  of  xi+X2+X3+xit=0, 
X3-hCi4+x5+xg=0,  Xg-fxg4x7+xg=0  and  X3+X3+X5+X7=0  over  the  finite 
field  GF(2)  excluding  the  treatments  (00  00  00  00)  and 
(11  11  11  11);  and  the  other  4  treatments  are  (10  00  00  00), 

(00  00  10  00),  (00  00  00  10),  (00  00  00  01). 

m=9 

The  treatment  matrix  consists  of  28  treatments  as  follows. 

The  treatments  are  solutions  of  xi+X2+X3+X4=0,  X3+X4-iX5+xg=0, 
xg4xg+X7+xg=0  and  X|+X3+X5+X7=0  over  the  finite  field  CF(2) 
excluding  the  treatments  (000000000),  (000000001), 

(11  11  11  11  1),  (11  11  11  11  0). 

m=10 

The  treatment  matrix  consists  of  35  treatments  as  follows. 

The  30  treatments  are  solutions  of  xi-ix2+x3+x4=^> 
X3+x4+x5-rxg=0,  X5+xg+x7+xg=0,  x 7+x g+x g+x | g  =0  and  x | +x 3+x 5+x7=0 
over  the  finite  field  CF(2)  exlcuding  the  treatments 
(00  00  00  00  00),  (11  11  11  11  11),  and  the  other  five  treatments 
are  (10  00  00  00  00),  (00  00  10  00  00),  (00  00  00  10  00), 

(00  00  00  00  10),  (00  00  00  01  00). 

The  arguments  for  the  cases  m=7 ,  8,  9,  10  are  similar  to  the 
argument  for  m=6. 
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4.  FINAL  REMARKS 

This  paper  touches  the  major  developments  in  factor  screening 
designs  over  the  period  of  more  than  25  years  and  deals  with  the 
important  issue  of  measuring  the  "dispersion"  effects  in  addition 
to  the  "location"  effects  of  factors  using  search  linear  models. 
The  influence  of  Professor  R.  C.  Bose  to  this  author  is  not  easy 
to  describe  in  words  and  this  author  takes  pride  in  dedicating 
this  work  to  honor  Professor  R.  C.  Bose. 
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